Draft version September 25, 2012 

Preprint typeset using L^-T^X style cmulateapj v. 5/2/11 



(N 

o 

(N 
D 

(N 



o 

u 

6 

(N 

> 
m 

00 

o 

(N 



13 



RECONSTRUCTING THE 7-RAY PHOTON OPTICAL DEPTH OF THE UNIVERSE TO z ~ 4 FROM 

MULTIWAVELENGTH GALAXY SURVEY DATA 

Kari Helgason 1,2 and Alexander Kashlinsky 2,3 

department of Astronomy, University of Maryland, College Park, MD 20742, USA; kari@astro.umd.edu 
2 Observational Cosmology Laboratory, Code 665, NASA Goddard Space Flight Center, Greenbelt MD 20771 and 
3 SSAI, Lanham, MD 20706; alexander.kashlinsky@nasa.gov 
Draft version September 25, 2012 

ABSTRACT 

We reconstruct 7-ray opacity of the universe out to z < 3 — 4 using an extensive library of 342 
observed galaxy luminosity function (LF) surveys extending to high redshifts. We cover the whole 
range from UV to mid-IR (0.15-25/xm) providing for the first time a robust e mpirical calcula t ion o f 
the 77 optical depth out to several TeV. Here, we use the same database as iHelgason et al.l (|2012T ) 
where the extragalactic background light was reconstructed from LFs out to 4.5/im and was shown to 
recover observed galaxy counts to high accuracy. We extend our earlier library of LFs to 25/iin such 
that it covers the energy range of pair production with 7-rays (1) in the entire Fermi/LAT energy 
range, and (2) at higher TeV energies probed by ground-based Cherenkov telescopes. In the absence 
of significant contributions to the cosmic diffuse background from unknown populations, such as the 
putative Population III era sources, the universe appears to be largely transparent to 7-rays at all 
Fermi/LAT energies out to z ~ 2 whereas becoming opaque to TeV photons already at z < 0.2 and 
reaching r ~ 10 at z = 1. Comparison of the currently available Fermi/LAT gamma-ray burst and 
blazar data shows that there is room for significant emissions originating in the first stars era. 
Subject headings: cosmic background radiation — galaxies: evolution — gamma rays: general 



1. INTRODUCTION 

The extragalactic background light (EBL) supplies 
opacity for propagating high energy GeV-TeV photons 
via an electron-positron pair production (77 — > e + e~) 
(jNikishov Ill962t ISteckerlll971I ). Determining the trans- 
parency of the universe is of fundamental importance 
for a wide variety of current observatories such as the 
space-borne Fermi/LAT instrument operating at ener- 
gies < 250 — 300 GeV to ground-based 7-ray telescopes 
probing energies >1 TeV. The distance at which the op- 
tical depth due to this interaction is r ~ 1 defines a 
horizon of the observable universe at 7-ray energies, and 
has been a subject of extensive efforts designed to model 
the build-up of EBL with time from the posited emis- 
sion h istory of galaxy populat i ons (e.g.. IStecker et all 
20061 iFranceschini et all l200l IKneiske fc Dole! 120101: 



Dommguez et al.ll2011l) 

In this Letter we show that, with observed galaxy 
populations over a wide range of wavelengths, one can 
uniquely reconstruct the optical depth of the universe at 
these energies out to redshifts z ~ 4. This empirical re- 
construction relies exclusively on data from an extensive 
library of galaxy luminosity functions (LFs) encompass- 
ing 18 finely sampled wavelengths from UV to mid-IR 
(0.15-24 /zm) relevant for the pair-production opacity. 
This methodology enables robust calculation of the 7- 
ray opacity in the Fermi/LAT energy range using galaxy 
surveys probing A < 4.5/im out to z < 4. Extending 
to TeV energies, probed by the ground-based Cherenkov 
observations, we use measurements out to 24 /xm; this 
extrapolation is robust for the redshifts currently probed 
these observations. This heuristic reconstruction using 
the observed galaxy populations defines the absolute floor 
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of the photon-photon optical depth due to known galaxy 
populations and deviations from it would allow the char- 
acterization of any emissions ina c cessible to dir ect tele- 
scopic studies (|Kashlinskvll2005bl: lGilmorell2012D. 

We use the methodology developed in IHelgason et all 
(|2012[) of reconstructing the EBL from observed galaxy 
populations in a compilation of 342 measured LFs cov- 
ering the UV, optical and near-IR; that compilation is 
slightl y updated compared to Table 1 of IHelgason et all 
(2012). The wealth of galaxy survey data has recently 
reached adequate redshift coverage to make such empir- 
ical estimation of the evolving EBL feasible and the re- 
construction was shown to reproduce independent data 
from g alaxy counts and t he cosmic infrared background 
(CIB, iKashlinskvl 12005a) . Our approach is completely 
independent of theoretical modeling describing the evo- 
lution of galaxy populations in that we use the LF data 
directly observed at all wavelengths out to z ~ 3 — 8 in 
this heuristic reconstruction from which we derive the op- 
tical depth due to pair production (see also IStecker et all 
|20H . 

Standard cosmological parameters are used below: 
Otot = 1, &m = 0.3, H = 70km • s^Mpc" 1 . 

2. RECONSTRUCTING THE EVOLVING EBL 
FROM DATA 

Quantifying the optical depth of the universe to high 
energy photons requires knowledge of the properties 
of the intervening EBL. Evolving galaxy populations 
compose the bulk of the EBL which is dominated by 
starlight in the UV/optical and thermally radiating 
dust at longer IR wavelengths. We use multiwave length 
survey data to fit the evolution of a single derived 
quantity, the luminosity density, in this otherwise 
assumption-free approach. We cover the whole range 
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Fig. 1. — Measurements of the evolving luminosity density in 
our binned rest-frame wavelength range 0.14-1.25/im in units of 
L©Mpc -3 . The solid curves are fits through the data points 
according to our fitting function (Equation r3{ with best-fit pa- 
rameters displayed in each panel as (log 10 a^,6^,c^) where 
is also in LqMpc~ 3 . All data points have been converted to 
units with h = 0.7. T he references f o r the data points can 
be found in Table 1 of Hclgason et al. (2012 ) with th e add i- 
tion of UV data from lSteidel et al.l 119991) ; ISullivan et alJ [120001); 
ITrever etall 120051); ISchiminovich et all (120051); IBudavari eTafl 
(2005); Sawicki & Thompson (2006); Finkclstcin et al. (20f2). 



from UV to mid-IR (0.15-25/im), providing for the first 
time an empiric al derivation of the 77 optical depth out 
to several TeV. iStecker et al.1 (|2012f) have used a similar 
approach to reconstruct the EBL at < 0.7/im. The 
survey data used in this Letter extends the collection of 
LFs presented in Hclga son et al.1 (120121) with expanded 
cover a ge in the UV (ISteidel et al.l 119991: ISullivan et al.l 



20001; ITrever et al.l (20051: ISchiminovich et all 



Budavari et al.1 2005; I S awicki fc Thompson! 
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Tresse et al J 120071: iFinkelstein et al.1 12012ft a nd out 
to mid-IR wave le ngths dRuiopakarn et al.l 120101: 
Shupe et all 119981: iXu et al.1 119981: ISanders et al l 
2003t ILe Floc'h et all 120051: IMagnelli et al.1 120111: 
Perez-Gonzalez et al. [|2005(). This library now contains 
342 measured LFs and allows us to reconstruct the 
evolving EBL and its spectrum in a finely sampled wave- 
length grid enc o mpass ing 0.15-24/mi out to z ~ 2 — 8. 
IHelgason et al.l (|2012D used this data to accurately 
recover the observed galaxy number counts in the 
0.45-4.5/im range. This therefore assumes only the 
existence of populations now observed out to z ~ 2 — 8 
and magnitudes as faint as toab ~ 23 — 26; additional 
populations, such as the hypothetical galaxies with first 
stars would then exist at still earlier times and have 
much fainter fluxes. 

In the rest-frame UV to near-IR, the LF is well de- 
scribed by the conventional Schechter parameterization 
whereas at mid-IR wavelengths, the LF seems to be 
better described by a broken power-law or a double- 
exponential profile. Regardless of the functional form, 
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Fig. 2. — Same as Figure [JJ for infrared wavelengths, 1.6- 
24//m. T he references fo r the data points can be found in 
Table 1 of Hcleason et al. (201 2]) with t h e add it ion of mid-IR 
data from [Shupe et all 11998): IXil et al.l (119981): ISanders et alJ 
(|2003|); ILe Floc'h et IT] 120051); IPerez-Gonzalez et alJ 120051 ); 
IRuiopakarn et al.l 1201(1) ; IMagnelli et al.l 1201*111 



the LFs can be integrated to give the comoving volume 
emissivity (we refer to this as the luminosity density) in 
the given rest-frame band 



e v (z) = / L v (j)(L l ,,z)dL v 



(1) 



where <j)(L Ul z) is the measured LF. Each data point in 
Figures Q] and [5] represents the luminosity density given 
by the autho rs along with la error bars (for references, 
see Table 1 in IHelgason et all (|2012l )). In the cases where 
this value is not given explicitly in the original papers, we 
have integrated the best-fit parameterized LF to obtain 
t v and have estimated the errors from the distribution 
of all of the values of e„ allowed within the la solutions 
of the individual fit parameters. However, mutual com- 
parison of uncertainties among the many different stud- 
ies is not very meaningful since some authors include 
various effects in addition to the statistical errors from 
the method of LF estimation, such as cosmic variance, 
fc-corrections, incompleteness and photometric system. 
Here, we have chosen to maximize our wavelength and 
redshift coverage by letting all available measurements 
contribute to our fitted evolution regardless of the er- 
ror treatment. In the cases where the median redshift 
of the sample is not explicitly given, we have placed the 
measurements at the midpoint of the redshift bin. Our 
wavelength interval is sampled at the rest-frame bands 
shown in the panels of Figure Q] and [2] where most of the 
LFs have been measured. The offsets from these defined 
wavelengths due to filter variations (e.g. SDSS vl and 

Johnson i7) are small enough to be negl ected. 

Mot ivated by the fitting formulae in [Helgason et al.1 
(|20lJ), which we found to be reliable over a wide range of 
wavelengths, we consider the following three parameter 
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Fig. 3. — Upper: the evolving EBL resulting from our empirical reconstruction. The different lines illustrate the build-up of EBL with 
cosmic time leading to the present day levels (thick solid curve). We display the evolving EBL in comoving coordinates i.e. without the 
expansion factor (1 + z) 3 for better mutual comparison. The left pointing triangles show the EBL threshold energy needed to interact 
with an observed 30GeV photon (blue), 300GeV, (red) and ITeV (green) originating at the redshifts shown. Lower: ou r reconstructed 
EBL c ompared to integrated counts in the literature along with the region bound by the upper/lower faint-end scenarios in Hclgason ct al. 
H2012I) (shaded). Th e black an d gray lines rep resent th e case of j? cuto ff = 10.2e V and 13.6eV respectively. The cou nts data are from 
^Gardner et aLfeoOCl) (purple) HXu et al.H2005l) (cyan) IIMadau fc Pozaettil 12000) (red), HKeenan et al.|[20lCT) (green). UFazio et alll200l) 
(orange ), jH opwood ct al. 2010) (cyan), HMetcalfe et al.ll2003j ) (black). l lPapovich et al.ll2004l ) (deep purple). I IBethermin et al.11201011 (light 
green), l lCharv et al.H2004l) (blue). 



fits for the evolution of the luminosity density 

e v (z) = a\(l + (z - z )) bx exp (-c\(z - z )) (2) 

where we fix zq = 0.8. Although not restricted to 
Schechter LFs, this functional form for e„(z) is never- 
theless equivalent to the underlying Schechter param- 
eters evolving as L* oc (1 + (z — z n )) bx and <f)* oc 
exp (— c\(z — zo)) with a constant faint-end slope. In- 
terpolating the rest-frame data between our 18 reference 
bands defines the rest-frame emissivity spectrum at any 
given epoch. We apply a cutoff to the spectrum above 
the Lyman limit, E cuto g — 13.6eV, corresponding to ef- 
ficient absorption of ionizing photons by hydrogen in the 
local environments. At all lower energies, the universe 
is assumed to be completely transparent to background 
photons. The integrated light from galaxies seen today 
is (e.g. Peebles 1971) 
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where v' = v{l + z) is the rest-frame frequency and 

dt 

dz 
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i/ (l + z)V(T 

In iHelgason etHI (|2012tl we considered two limiting 
cases for the evolution of the faint-end of the LF and 
showed that the distribution of galaxies from LF data, 



when projected onto the sky, accurately recovered the 
observed galaxy counts across the optical and near- 
infrared. The flux from the integrated counts is dis- 
played as shaded regions in the lower panel of Figure [3] 
along with our empirically determined EBL (solid line) 
which is in good agreement with integrated counts data 
in the literature apart from wavelengths > 6/im where 
the steep evolution of the mid-IR LFs (8-24/xm) causes 
our EBL to be a factor o f ~2- 3 higher than integrate d 
counts fromlCharv et al.1 (12004D . iHopwood et all (igOlOD . 
iPapovich et all (|2004f ) and iBethermin et all (|2010D . Al- 
though these authors do not claim to fully resolve 
the CIB at these wavelengths, the discrepancy is large 
enough to indicate a mismatch between number counts 
and mid-IR LF measurements at z > 0.5. This issue is 
apparently also encoun tered in galaxy evolution models; 
ISomerville et al.l (|2012f) are not able to simultaneously 
account for the integrated counts and the bright-end of 
the observed LF in the 8-24/im range using different dust 
templates. In fact, this is the wavelength regime where 
varying degrees of dust contribution and PAH emission 
make the spectrum less predictable. Recent upper limits 
derived from a TeV source s pectra also favo r low levels of 
CIB at these wavelengths (|Orr et all 1201 It iMever et al.l 
2012). At this stage, one must therefore question the 
robustness of EBL reconstructed from LFs at > 8/im. 

The UV/blue end of the EBL turns out to be sensitive 
to the abundance of photons with energies just below 
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the Lyman limit, 13.6eV. The redshifted far-UV contri- 
bution dominates the EBL below 0.5/xm due to the steep 
increase of the star formation rate at earlier times. For all 
galaxy types, there is considerable absorption in the Ly- 
man series which we do not account for and we illustrate 
this dependence by considering the case where Lyman- 
series absorption completely suppresses the spectrum 
above 10.2eV (instead of E cuto s — 13.6eV); shown as 
gray lines in Figure[3] (lower panel). We subsequently dis- 
play our optical depths for the both cases which bracket 
the true behavior. 

3. THE PHOTON-PHOTON OPTICAL DEPTH 

The relevant quantity for computing the optical depth 
due to photon-photon interaction is the rest-frame num- 
ber density of photons as a function of time and energy, 
n(E, z). We shall refer to the energy of a photon belong- 
ing to the EBL as E and we use £ for the propagating 7- 
rays. Rest-frame quantities are denoted with a prime. At 
any given epoch, the photon number density (in proper 
coordinates) is composed both of sources emitting in the 
rest-frame as well as the accumulated emission from ear- 
lier times 



n(E',z) = (1 + zf 



e u ,(z')/h dt 
hu 1 ~dz~' 



dz' 



(5) 



where h is the Planck constant (the extra h is to con- 



vert e v to per unit energy, dE — hdv) and 



(1 



z')/(l + z). The condition for pair production is that 
the total energy in the center-of-mass frame must sat- 
isfy £'E'(1 — cosO) > 2(m e c 2 ) 2 where 9 is the angle of 
incidence. This means that in order to interact with a 
7-ray of energy £', background photons must have wave- 
lengths of < 1.0(5' /210GeV)^m. The cross section for 
this interaction is 

a{E',£',^) = ^{l-^) 



2p{(3 2 -2) + (3-/? 4 )ln 



where 



1-/3 
(6) 



T- 



2m 2 c 4 



E'£'(l-fi)' 



[JL = COS t 



For the most likely angle of incidence, fi = (side-on), 
the probability for interaction is maximized at roughly 
four times the minimum threshold energy, ~ 4m e c 2 /£. 
The optical depth encountered by a high energy photon 
originating at z (r ~ anl) can be expressed in terms of 
its observed energy, £, as 



J£,z) 



dt 



a 



d^i 



: / a(E',£',fj,)n(E',z')dE' 

/2m 2 c 4 /£'(l-^) 

(7) 

where n(E, z) comes from Equation [5] In Figure |4] we 
display the calculated optical depths as a function of ob- 
served energy for 7-rays originating at several redshifts. 
The optical depth roughly traces the shape of the num- 
ber density of EBL photons with a sharp drop in the 
optical depth at the lowest energies due to the cutoff at 
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Fig. 4. — Solid lines show the 77 optical depth contributed 
by known galaxy populations assuming -E cu t ff = 13.6eV (gray 
curve in the lower panel of Figure [3j . The curves are not drawn 
beyond the energy of (m e c 2 ) 2 / i?24^m(l + 2) as we do not con- 
sider data at A > 25/xm. The shaded regions show the bound- 
ar ies of the upper / lower scenarios in the empirical reconstruction 
of IHelga son ct al. (2012) out to 4.5/xm. The dotted vertical line 
shows roughly the highest energy probed by Fermi/LAT and the 
dashed line shows where r = 1 for reference. 

the Lyman limit. We also show the regions encompassed 
by the two scena rios of the faint-end ev olution from the 
reconstruction of Hclg ason et al.l ()2012D . 

4. APPLICATION TO HIGH ENERGY 
OBSERVATIONS 

What do these reconstructed opacities imply for obser- 
vations of high-energy sources with the current instru- 
ments? Blazars and gamma-ray bursts (GRBs) are ex- 
" imples of high energy extragalactic sources whose spec- 
tra is affected by attenuation of photons in excess of 
> lOGeV. Two types of datasets are relevant for this 
discussion: space-borne Fermi/LAT measurements at 
< 300GeV and ground-based telescopic measurements 
extending out to ^TeV energies. Because extragalactic 
7-ray absorption increases with both with redshift and 
energy, the EBL can be constrained based on the highest 
energy photon s observed from a source provided the red- 
shift is known (jAbdo et al. 2010). GRBs have the advan- 
tages of being observable across great distances and typi- 
cally displaying harder spectra than most blazars at sub- 
GeV energies. Figure [3] shows that for Fermi-observed 
sources it is sufficient to use data out to ~ 4.5/mi (red 
symbols), whereas for TeV range observations, survey 
data are needed out to the longer wavelengths (green 
symbols). Figure|3]shows the reconstructed optical depth 
explicitly confirming this. We now briefly discuss the im- 
plications: 

Fermi/LAT detects blazars and GRBs out to ener- 
gies ~ 250 — 300 Gev. In Figure [S] (upper panel) we 
show curves of constant 77 optical depth in £ — z space 
and compare with the most constraining high energy 
Fermi/LAT sourc es with known redshifts (taken from 
lAbdo et al.l l2010f) . The contours of r = 1,2,3,5 cor- 
respond to probabilities of photon being absorbed by 
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Fig. 5. — Upper: the curves show where a photon of energy £ 
originating at z, encounters exactly t 77 = 1,2,3,5 (solid, dotted, 
dashed, dash-dotted respectively). Black and gray correspond to 
the cases of the rest-frame spectrum cutoff energy is at 10.2 and 
13.6eV respectively. The symbols show the highest energy photon 
observed in GRBs (blue squares) and a selectio n of the most con- 
straining of Fermi/LAT blazars (orange circles) (Abdo ct aI1 l20101 ; 
IMcConviUe et aXll2011ll . Lower: the observed energy spectrum of 
the sources 1ES 1101-232 (z = 0.186; left) and 3C 279 (z = 0.536; 
right) shown with black diamonds, and the corresponding deab- 
sorbed data as blue squares. The deabsorbed spectra have best-fit 
photon indices r in = 1.49 and r- ln = 2.28 respectively, but deviate 
substantially from a power-law at the highest energy bins. 



z = 0.186 interacts mos t strongly with optical and near- 
IR background photons ([Aharonian et al. 2006]) . The ob- 
served spectrum is relatively hard (T = 2.88 ± 0.17) and 
results in a best-fit intrinsic photon index of r; nt = 1.49 
after deabsorption. For this particular source, the upturn 
at TeV energies is largely driven by the EBL photons at 
> 5/im and would be less pronounced if our EBL re- 
flected the integrated counts data in the lower panel of 
Figure [3l Because the EBL changes with time in both 
shape and amplitude, the effects on 7-ray absorption 
become even more prominent for more distant sources. 
The spectrum of the distant radio quasar 3C 279 has 
been captured by MAGIC during d i fferent flaring event s 
(jMAGIC CoUab'orationiraIll200l lAleksic et all 120 111) . 
The deabsorbed spectrum of 3C 279 shown in Figure [5] 
(lower) also deviates substantially from a simple power- 
law in the highest energy bins which is unlikely to be due 
to our near-IR background being overestimated. This 
behavior of th e deabsorbed 3C 279 sp ectrum has been 
pointed out by iDommguez et ail (|2011l) who suggest ei- 
ther improved emission models or instrumental system- 
atic uncertainties as potential solutions. Another possi- 
bility is that some fundamental effects are missing such as 
secondary 7-rays produced along the li ne of sight by cos- 
mic r ays accelerated by the blazar jet ()Essev fc Kusenkol 
2010). In absence of such secondary effects however, the 
universe should be completely opaque for TeV sources at 
z>0.5. 

Upper limits for the EBL derived from TeV spec- 
tra rely on assum ptions of the hardne ss of the intrinsic 
blazar spectrum. iMever et al.l (|2012D derive limits for 
the whole range of optical to far-IR EBL using an exten- 
sive source sample from both Fermi/LAT and ground- 
based Cherenkov telescopes. Their results allow a total 
integrated NIR flux (1-10/xm) of ~20nWm _2 sr _1 in ex- 
cess of known galaxy populations whereas there is, at 
most, little room for extra contribution in the mid-IR 
(> 10/xm). In fact, our LF-deri ved EBL i s incon sistent 
with the lowest mid-IR limits of lOrr et all (|2011l ). 



the EBL of 63%, 86%, 95%, 99.3%, respectively. In the 
absence of new populations, the universe remains fairly 
transparent at the Fermi/LAT energies out to z ~ 2 — 3. 
Our reconstructed EBL is fully consistent with all the 
available LAT data and, in fact, allows for non-negligible 
extra levels of the CIB from new populations such as pos- 
sibly have existed at higher z. As the Fermi mission pro- 
gresses and Figure[S] (upper) becomes more populated at 
the highest energies, sources at high-z will provide better 
constraints for the optical/NIR EBL. 

Ground-based Cherenkov telescopes have pro- 
duced good quality spectra for TeV-blazars, although for 
sources at significantly lower redshifts than Fermi/LAT. 
If the evolving EBL is known to a good accuracy, one 
can deabsorb observed blazar spectra to reveal the in- 
trinsic component, which is expected to have a power- 
law form, dN/dE oc E~ T . The lower panels in Figured] 
demonstrate how our reconstructed EBL affects the spec- 
trum of two known blazars at relatively high redshifts, 
both of which have been used to place upper limits on 
the optical/NIR EBL. Good quality spectrum of the BL 
Lac object 1ES 1101-232 (z = 0.186) has been obtained 
by HESS in the energy range 0.16 — 3.3TeV which at 



5. SUMMARY 

We have shown that it is possible to robustly recon- 
struct the evolving EBL in the universe using our earlier 
library of multiwavelength survey data now updated to 
exten d from the UV out to the mid-IR (Helgason et al. 
2012). This reconstruction uniquely defines the 7-ray 
opacity out to TeV energies for sources at z < 4 and 
shows that at the energy bands probed by Fermi/LAT, 
the universe is fairly transparent out to z ~ 2 — 3, un- 
less unknown sources at high redshifts contribute non- 
negligible amounts of CIB. At TeV energies, probed by 
ground based telescopes, the universe becomes optically 
thick at z ~ 0.5 so any such photons associated with the 
sources at higher redshifts would have to be of secondary 
origin. 

Our reconstructed EBL and optical depths are avail- 
able upon request. This work was supported by NASA 
Headquarters under the NASA Earth and Space Sciences 
Fellowship Program - Grant NNX11AO05H. KH is also 
grateful to The Leifur Eiriksson Foundation for its sup- 
port. We thank W. McConville, B. Magnelli and M. 
Ricotti for useful communications. 
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